IDA 


INSTITUTE  FOR  DEFENSE  ANALYSES 


Double  Negative  Materials  (DNM), 
Phenomena  and  Applications 


John  Franklin 
John  Biddle 
Bohdan  Balko 


July  2009 

Approved  for  public  release; 
distribution  is  unlimited. 

IDA  Document  D-3887 
Log:  H  09-000915 


This  work  was  conducted  under  IDA’s  central  research  program, 
CRP  2117.  The  publication  of  this  IDA  document  does  not  indicate 
endorsement  by  the  Department  of  Defense,  nor  should  the  contents  be 
construed  as  reflecting  the  official  position  of  that  Agency. 

©  2009  Institute  for  Defense  Anaiyses,  4850  Mark  Center  Drive, 
Alexandria,  Virginia  22311-1882  •  (703)845-2000. 

This  material  may  be  reproduced  by  or  for  the  U.S.  Government. 


INSTITUTE  FOR  DEFENSE  ANALYSES 


IDA  Document  D-3887 


Double  Negative  Materials  (DNM), 
Phenomena  and  Applications 


John  Franklin 
John  Biddle 
Bohdan  Balko 


Preface 


This  work  was  supported  by  a  Central  Researeh  Program  (CRP)  task  at  the  Insti¬ 
tute  for  Defense  Analyses  (IDA). 


Foreword 


When  we  started  the  work  on  this  doeument  in  Fiseal  Year  (FY)  2006,  the  double 
negative  materials  (DNMs)  or  negative  index  of  refraction  (NIR)  materials  research  com¬ 
munity  was  involved  in  justifying  basic  physical  questions  of  causality,  validity  of  the 
concept  of  negative  index  of  refraction  interpretation  of  experimental  results,  and  so 
forth.  At  the  same  time,  progress  was  being  made  in  developing  different  ways  of  con¬ 
structing  DNMs,  increasing  the  range  of  frequencies  toward  the  optical/visible  regime, 
approaching  three-dimensional  (3D)  operability,  reducing  attenuation,  and  widening  the 
frequency  window.  However,  even  now  in  2009,  questions  about  the  validity  of  the  DNM 
concept  are  still  being  challenged.  ^  All  this  has  been  going  on  while  practical  applications 
are  moving  forward  and  new  designs  of  metamaterial  structures  are  being  introduced  and 
tested. 

As  a  result  of  this  study,  new  areas  of  interest  to  IDA  opened:  (1)  superlenses  or 
beating  the  diffraction  limit,  (2)  reverse  Casimir  effect  for  reducing  friction/stiction  in 
microelectromechanical  systems  (MEMS),  (3)  cloaking  in  two  and  three  dimensions,  and 
(4)  the  development  of  transformation  optics.  We  believe  that  these  new  areas  of  interest 
are  worthy  of  investigation  by  the  Department  of  Defense  (DoD). 


1  Munk,  B.  A.  Metamaterials:  Critique  and  Alternatives  (New  York:  John  Wiley  and  Sons,  2009). 
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Executive  Summary 


In  this  report,  we  will  diseuss  the  nature  of  negative  index  of  refraetion  (NIR) 
materials  or  double  negative  materials  (DNMs),  their  realization  (i.e.,  how  they  are  pro- 
dueed),  their  speeial  effeets  and  potential  applieations,  and  their  modeling — starting  with 
simple  basie  models  and  eulminating  with  eomplex  eomputer  eodes. 

“Normal”  optieal  materials  that  have  eleetrie  permeability  £  >  0  and  magnetie 
permittivity  /r  >  0  are  quite  familiar  to  students  of  physios  and  engineering.  DNMs  that 
have  eleetrie  permeability  £  <  0  and  magnetie  permittivity  //  <  0  are  quite  different.  They 
do  not  normally  oocur  in  nature,  are  man-made,  and  have  some  unusual  properties.  The 
purpose  of  this  paper  is  to  explore  this  seoond  possibility. 

DNM  or  left-handed  material  (LHM),  although  the  subjeot  of  muoh  speoulation  in 
1964,  was  not  fabricated  until  the  year  2000  (at  a  narrow  bandwidth  centered  around 
4.84  GHz).  The  term  “metamaterial”  was  used  by  the  Defense  Advanced  Research 
Projects  Agency  (DARPA)^  for  these  new  materials.  These  materials  are  constructed 
from  normal  materials  (e.g.,  copper)  in  such  a  way  that  they  act  as  scattering  centers  with 
special  properties  for  affecting  electric  and  magnetic  fields.  The  novel  bulk  electromag¬ 
netic  (EM)  properties  emerge  from  the  collective  response  of  subwavelength  structures 
made  out  of  the  normal  structures.  Since  the  year  2000,  many  other  DNM  have  been  con¬ 
structed  using  various  techniques  and  pushing  the  operational  region  into  higher  frequen¬ 
cies,  including  the  visible. 


2  In  2001,  the  Defense  Advanced  Research  Projects  Agency/Defense  Sciences  Office  (DARPA/DSO) 
set  up  a  metamaterials  research  program.  DARPA’s  initial  interest  was  in  applications  beyond  optics, 
such  as  magnets  for  more  powerful  motors,  radar  applications,  and  textured  surfaces  that  could  control 
a  material’s  heat-handling  capabilities.  DARPA  later  extended  its  interest  to  advances  in  ultra-high- 
frequency  antennas  for  possible  radar  applications. 
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1.  Introduction 


1,1  Double  Negative  Material  (DNM) — ^What  Is  It? 

Recently,  much  interest  has  been  focused  on  metamaterials  that  show  unusual 
interactions  with  electromagnetic  (EM)  waves — interactions  that  have  not  been  seen 
before  with  normal,  naturally  occurring  materials.  These  materials  have  been  referred  to 
as  negative  group  velocity  (NGV)  materials,  negative  index  materials  (NIMs),  double 
negative  materials  (DNMs),  left-handed  materials  (LHMs),  and  negative  index  of  refrac¬ 
tion  (NIR)  materials.  All  these  names  are  appropriate,  but,  in  this  paper,  DNM  will  be 
used  most  often. 

EM  waves  in  space  contain  an  electric  field  (E)  and  a  magnetic  field  (H)  perpen¬ 
dicular  to  the  direction  of  propagation,  k,  as  shown  in  Eigure  1-1. 


Figure  1-1.  EM  Fields  in  a  Plane  Propagating  Wave 

Propagation  of  EM  waves  is  described  in  terms  of  the  phase  velocity,  given  by 
vp  =  — ,  where  X  is  the  wavelength  and  T  is  the  period  of  the  wave,  and  the  group  veloc- 

ity,  given  by  vg  =  — ,  where  co  is  the  frequency  of  the  wave  and  k  is  its  wave  number.  3 
8k 

The  energy  flow  in  space  is  given  by  Poynting’s  vector,  S  =  E  x  H.  The  two  velocities 
describe  two  different  properties  of  a  wave:  the  phase  velocity,  Vp  ,  with  which  the  phase 
of  a  wave  of  one  frequency  travels  and  the  group  velocity,  Vg  ,  with  which  the  envelope 


3  Sometimes  the  parameters  Vp,  Vg,  and  k  will  be  bolded  to  indicate  that  they  are  vectors  and  have  a 
direction  associated  with  them. 
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of  the  wave  propagates  and,  in  most  cases,  can  be  thought  of  as  the  velocity  of  energy 
transfer. 


Figure  1-2  shows  the  difference  between  group  and  phase  velocity  in  normal 
materials  (i.e.,  right  handed  materials  (RHMs))  and  DNMs  or  LHMs). 


Group  vs.  Phase  Velocity 

For£>0,  |j>0  =^(k,  E,  H)  right-handed  set 


k  s 


''gTTVp 


For£<0,  |j<0  =>  (k,  E,  H)  left-handed  set 


H 


Figure  1-2.  Orientation  of  Vectors  in  Positive  index  of  Refraction  (PIR)  Materiai  and  DNM 
Note  1  for  Figure  1  -2:  /n  the  figure,  u  is  the  energy  density,  Vg  is  the  group  veiocity,  and  Vp  is 
the  phase  veiocity.  Vg  is  paraiiei  to  S  ,  and  Vp  is  paraiiei  to  k  . 

Note  2  for  Figure  1-2:  kxE  =  O)  juH ,  and  S  =  ExH  =  {u)vg  . 


The  simplest  example  of  the  difference  between  normal  materials  and  DNMs  is 
obtained  by  comparing  the  effects  of  the  two  materials  on  refraction.  In  a  normal  trans¬ 
parent  (dielectric)  material,  an  EM  beam  is  refracted  toward  the  normal  (blue  line)  as  it 
passes  from  air  into  the  higher  optical  density  (positive  index)  material  as  shown  in  Fig¬ 
ure  1-3.  In  a  DNM,  the  beam  is  refracted  to  the  other  side  (red  line)  of  the  normal  as 
shown  in  Figure  1-3. 

The  angle  of  refraction  is  obtained  from  Snell’s  law,  which  states  that  the  ratio  of 
the  sines  of  the  angles  of  incidence  and  refraction  is  equal  to  the  ratio  of  velocities  in  the 
two  media  or,  equivalently,  to  the  inverse  ratio  of  the  indexes  of  refraction  ni  and  n2. 
Thus, 

nisin(6i)  =  n2sin(62),  (1-1) 

where  6i  and  62  are  the  angles  of  incidence  and  refraction  respectively,  as  shown  in  Fig¬ 
ure  1-3.  For  a  DNM,  the  law  still  holds.  Since  112’  <  0,  and  62’  <  0,  the  product  n2 ’sin(62 ) 
is  positive,  but  the  refracted  beam  direction  is  to  the  left  since  62’  <  0,  as  shown  in  Fig¬ 
ure  1-3. 
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Figure  1-3.  Refraction  in  Normal  Material  (Blue  Line)  and  DNM  (Red  Line) 

Figure  1-4  shows  the  practical  application  of  this  DNM.  In  this  figure,  we  com¬ 
pare  the  effects  of  focusing  the  light  of  lenses  made  from  normal  (blue)  and  DNM  (red). 
Notice  that  the  DNM  has  the  opposite  effect  from  that  of  the  normal  material:  convex 
lenses  diverge  the  beam  while  concave  lenses  converge  the  beam.  In  Figure  1-5,  a  plate 
containing  normal  material  (a)  diverges  the  beam  (si)  so  that  a  virtual  source  is  seen 
inside  the  material  and  a  plate  containing  DNM  (b)  focuses  the  beam  (si)  inside  the 
material  (fl)  and  on  the  other  side  from  the  source  outside  the  material  (fZ). 


(a)  (b) 

Figure  1-4.  Effect  on  EM  Beams  of  Lenses  Made  From  (a)  Normal  Material  and  (b)  DNM 
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Figure  1-5.  Effect  on  EM  Beams  of  Plane 
Plates  Made  From  (a)  Normal  Material  and  (b)  DNM 

1.2  Early  History 

Most  authors  attribute  the  first  theoretieal  study  of  DNM  to  V.  G.  Veselago 
(1967),  who  developed  the  basic  theoretical  underpinnings  of  the  interaction  of  EM 
waves  with  these  materials  and  discussed  various  unusual  effects  on  well-known  pheno¬ 
mena  such  as  refraction,  Snell’s  law  of  refraction,  the  Doppler  Effect,  Cherenkov 
radiation,  and  so  forth.  Others  preceded  him  in  considering  the  idea — if  not  for  electro¬ 
magnetism,  then  for  acoustics:  D.  V.  Sivukhin  (1957);  G.D.  Malyuzhinets  (1957);  E.I. 
Mandel’shtam  (1950);  H.  C.  Pocklington  (1905);  and  H.  Eamb  (1905),  who  may  have 
been  the  first  to  suggest  backward  waves  in  mechanical  systems.  However,  Veselago 
provided  a  theoretical  paper  that  had  the  first  serious  and  complete  discussion  of  DNM.^ 

1.3  Recent  History 

Whereas  the  early  history  concentrated  on  the  theoretical  analysis  of  various 
effects  within  the  realm  of  physical  possibility,  recent  history  has  been  concerned  with 
experimental  developments:  construction  of  DNMs  (metamaterials),  testing,  and  the 
search  for  useful  applications.  This  new  focus  does  not  mean  that  theoretical  studies  have 
been  abandoned — only  that  the  focus  is  on  applications  since  the  validity  of  the  pheno¬ 
mena  has  been  well  established.  The  experimental  verification  of  the  phenomena  by  R.  A. 
Shelby,  D.  R.  Smith,  and  S.  Schultz  in  2001  and  Andrew  A.  Houck,  J.  B.  Brock,  and  1.  E. 


^  Veselago,  V.  G.  “The  Eleetrodynamics  of  Substances  With  Simultaneously  Negative  Values  of  s 
and  p,.”  Soviet  Physics  Uspekhi  10  (4)  (January-February  1968):  509-514.  Available  at 
http://www.iop.Org/EJ/article/0038-5670/10/4/R04/PFlU  10  4  R04.pdf?request-id=be94154d-cd29- 

4796-aca6-209e7c9b79bb 
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Chuang  in  2003  established  this  trend.  The  hrst  experimental  evidenee  for  DNM  is 
shown  in  Figure  1-6.  The  experiment  was  performed  at  100  MHz,  with  two  prisms  (one 
made  of  Teflon  (normal  material)  and  the  other  made  of  a  DNM  (i.e.,  LHM)).  The  Teflon 
beam  is  refracted  at  30  deg  while  the  LHM  material  refracts  at  a  negative  angle  (in  this 
case,  -60  deg),  as  was  expected  with  a  DNM. 


Figure  1-6.  First  Experimentai  Evidence  for  DNM 

Source:  Shelby,  R.  A.,  D.  R.  Smith,  and  S.  Schultz.  April  6,  2001.  “Experimental  verification  of  a 
Negative  Index  of  Refraction.”  Science  292  (5514):  77-79. 

Notel  for  Figure  1-6:  DNM  result  solid  (blue  line)  line;  normal  material  dashed  line  (green  line). 
Note  2  for  Figure  1  -6:  The  LHM  consist  of  square  copper  split  ring  resonators  and  copper  wire 
strips  on  fiberglass  circuit  board  material.  The  rings  and  wires  are  on  opposite  sides  of  the 
boards,  and  the  boards  have  been  cut  and  assembled  into  an  interlocking  lattice.  The  LHM  had  a 
unit  cell  dimension  of  5  mm,  a  factor  of  6  smaller  than  the  X-band  (8  to  12  GHz)  center  wave¬ 
length  of  3  cm. 

1.4  Trends 

The  chart  in  Figure  1-7  shows  the  number  of  citations  in  the  literature  to  Vese- 
lago’s  seminal  paper  on  DNM  from  2000  to  2003.  In  Figure  1-8,  we  show  the  number  of 
publications  on  LHMs  from  1990  to  2006  obtained  by  a  “Google  Scholar”  search. 

1.5  Supplementary  Material 

Appendix  A,  which  is  on  a  compact  disk  (CD)  attached  to  the  inside  back  cover  of 
this  document,  is  a  briefing  on  DNMs  that  was  presented  in  October  2008.  This  briefing 
contains  movie  clips  (pages  A-7,  A-29,  and  A-30)  that  show  a  simulation  of  EM  waves 
propagating  through  normal  material  and  through  DNM. 
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Figure  1-7.  Number  of  Citations  of  Veselago’s  Seminal  Paper  on  DNM 

Source:  MITRE  Corporation.  November  2003.  JASON  Report  on  Metamaterials.  MITRE 
Technical  Report  03-510.  McLean,  VA.  MITRE  Corporation. 


Figure  1-8.  Papers  on  DNM  From  1990  to  2006 
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2.  Background 


To  understand  the  DNM  phenomena  and  the  exeitement  it  is  eausing  and  to 
appreeiate  the  many  new  applieations  that  researehers  are  eonsidering,  we  need  to 
examine  the  response  funetions  of  materials  to  EM  fields.  To  help  supplement  this  inves¬ 
tigation,  Appendix  B  presents  a  glossary  of  terms  used  in  DNM  diseussions,  and  Appen¬ 
dix  C  eontains  a  bibliography  of  DNM. 

Two  maeroseopie  parameters  that  deseribe  materials  and  their  interaetion  with 
EM  fields  are  the  eleetrie  permittivity,  s,  and  the  magnetie  permeability,  ju.  The  parameter 
s  determines  how  well  a  material  polarizes — or  ereates  an  eleetrie  dipole,  P.  The  permit¬ 
tivity  defines  the  displaeement  D  =  sE  =  sqE  +  P.  The  parameter  ji  determines  how  well 
a  material  magnetizes — or  ereates  a  magnetie  dipole,  M — and  defines  the  magnetization 
H  =  B/ji  or  B  =  jiH  =  jUoH  +  M. 

The  index  of  refraetion,  n,  is  usually  ealeulated  as  the  produet  (jie)  ).  Eor  nega¬ 
tive  ji  and  e,  a  simple  replaeement  of  ji  with  -//  and  e  with  -e  still  gives  a  positive  value 
for  the  index  of  refraetion  sinee  n  =  ((-ji  )(-e  )/  =  (jie)  is  still  positive.  However,  the 
square  root  has  two  solutions:  one  positive  and  one  negative.  Thus,  the  simple  solution 
when  fi  and  s  are  both  negative,  is  to  take  the  negative  square  root  to  get  the  index  of 
refraetion: 

n  =  (±)  (so  pick  n=  (2.1) 

However,  we  have  a  more  satisfaetory  solution  in  eomplex  spaee.  Elsing  eomplex 
notation,  the  result  eomes  out  more  naturally  with  -p  expressed  as  p  and  -s  expressed 
as  E  so  that  n  =  [ (-p)  (-e)]‘^^  (p  e  =  (p  =  -(pe/'^^. 

An  EM  wave  propagating  through  a  medium  is  eharaeterized  by  its  wave  veloeity 
and  group  veloeity  and  its  interaetion  with  the  material  through  material  properties  p  and 
£.  The  relationship  between  the  wave  veetor  k  and  the  veetors  E  and  H  is  given  by 

kxE  =  CO  juH  (2.2) 

for  a  sinusoidal  wave  of  frequeney  co. 
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The  Poynting  vector,  S,  represents  energy  flow  and  is  given  by 
S  =  E  X  H  =  —E  X  B  ={u)  Vg 

where  u  is  the  energy  per  volume  and  the  group  velocity. 

Figure  1-1  shows  the  difference  between  group  and  phase  velocity  in  normal 
materials  (i.e.,  RHMs)  and  DNM  or  LHMs.  How  an  EM  beam  propagates  through  a 
material  and  how  it  passes  through  an  interface  between  two  materials  of  different  optical 
densities  are  determined  by  these  two  parameters  (i.e.,  group  and  phase  velocity).  We 
examine  this  question  next  using  Maxwell’s  equations. 

2.1  Maxwell’s  Equations 

To  study  materials  and  their  interaction  with  EM  fields  theoretically,  a  useful 
starting  point  is  Maxwell’s  equations,  which  govern  the  dynamics  of  EM  fields  and 
charges.  In  Maxwell’s  equations,  both  //  and  e  appear,  as  does  the  index  of  refraction, 
which  describes  the  optical  density  of  a  material  and  appears  as  a  combination  of  the  two 
in  the  form  n  =  (jis)  .  In  the  usual  treatments,  we  assume  that  the  result  is  a  positive 
number;  however,  it  has  been  determined  recently  that  this  result  (a  positive  number)  is 
not  the  only  possibility  (see  Eqs.  (2.4)-(2.7)). 

In  the  meter,  kilogram,  and  second  (mks)  system  of  units  and  in  a  region  without 
free  charge.  Maxwell’s  equations  are  as  follows: 


y-£E  =  o, 

(2.4) 

V-juH  =  0, 

(2.5) 

dt 

(2.6) 

and 

V7  TT 

VxH  =  —  ,  (2.7) 

ot 

where  E  is  the  electric  field  intensity  in  units  of  volts  per  meter  and  H  is  the  magnetic 
field  in  units  of  amperes  per  meter.  In  these  units,  s  =  -8.85  pE/m  and  the  permeabil¬ 

ity  is  /r  =  /To  =  1 .25  pHz/m.  The  equations  could  also  be  written  using  D,  the  electric  flux 
density  in  units  of  coulomb  per  square  meter,  and  B,  the  magnetic  flux  density,  in  units  of 
Weber  per  square  meter. 
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Assuming  isotropic  e  and  /r,  we  can  obtain  the  wave  equation  in  the  two  forms: 

W^E  =  siu^  (2.8) 

dr 

and 

=  (2.9) 

dr 

-1/7 

Clearly,  (jue)~  is  the  wave  velocity  in  both  forms.  The  flow  of  energy  is  given  by 

S  =  E  X  /T,  where  S  is  the  Poynting’s  vector  in  units  of  Joules  per  meter  -sec  and  the 

2  2 

energy  density  is  m=  (1/2)  (eE  +  jiH  ). 

For  demonstrating  some  basic  phenomena,  let  us  assume  plane  waves — specifi¬ 
cally,  a  linear  polarized  beam  with  £■  in  x  direction  {EL  in  y  direction),  moving  in  the  z 
direction  with  field  sources 


E 


..  i{k-z-cot) 

2C'  ^ 


(2.10) 


and 


H  = 


(OjLl 


i{k-z-CL)t)  + 

■y-e 


n 

2 
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(2.11) 
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where  k  =  k^  and  k  =  ejua) .  The  fields  are  perpendicular  to  each  other  and  in  a  plane 
normal  to  the  direction  of  propagation  k  as  shown  in  Figure  1-1  in  a  vacuum.  We  have 
introduced  a  phase  difference  of  nil  between  the  E  field  (e.g.,  Eq.  (2.10))  and  the  H  field 
(e.g.,  Eq.  (2.11)),  realizing  that  although  Eq.  (2.8)  and  Eq.  (2.9)  are  decoupled,  the  fields 
are  actually  related  energy  fields.  The  phase  difference  between  the  two  fields  makes  the 
energy  of  the  photon  represented  here  constant  during  the  whole  cycle.  The  electric  field 
feeds  the  magnetic  field  and  vice  versa. 


Since  we  are  only  interested  in  real  parts  of  E  and  H,  evaluating  the  time  averaged 
Poynting  vector  S  and  the  energy  density  U  gives 

1 

{S)= - -k  (2.12) 

coin  2 


and 
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(2.13) 


_  {SE^)  ^  _ 
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From  Eq.  (2.12),  we  see  that  when  ju  >  0,  S  and  k  (energy  and  phase)  are  in  the 
same  direetion),  but,  when  ju  <  0,  S  and  k  (energy  and  phase)  are  in  the  opposite  diree- 
tion.  This  differenee  is  a  major  distinetion  between  normal  materials  and  DNM  and  is  the 
basis  for  later  diseussion. 

Eq.  2-13  for  the  energy  density  applies  to  vaeuum  and  PIR  materials.  Eor  disper¬ 
sive  materials,  the  equation  derived  by  Eandau  and  Lifshitz^  and  used  by  Veselago®  is 


dco  dco 


(2.14) 


Eor  the  energy  density  to  be  positive,  the  requirement  on  the  derivatives  is 


S{sa))  ^  Q  dipiQ})  ^  Q 


(2.15) 


dco  ’  dco 


These  requirements  allow  e  and  ji  to  be  simultaneously  negative  if  they  are  functions  of  <3? 
(i.e.,  e(a)  and  ju(co)). 

2.2  Possible  Combinations  of  e  and  in  Materials 

Now,  if  we  consider  Eq.  (2.8)  and  Eq.  (2.9)  and  look  at  all  the  possible  combina¬ 
tions  of  E  and  ju,  we  find  that  a  propagating  wave  is  supported  by  the  material  only  when 
jU  >  0,  £  >  0  and  when  ji  <  0,  e  <  0.  When  jj.  >  0,  e  <  0  or  when  ju  <  0,  e  >  0,  the 
propagating  wave  decays  and  is  not  supported  by  the  material,  as  shown  in  the  /j.-e  plot  in 
Eigure  2-1.  Since  no  propagation  can  occur  in  these  regions,  they  are  not  of  interest  here. 
When  the  product  e/u  is  positive,  we  get  transmission  across  a  boundary.  When  the  prod¬ 
uct  qu  is  negative,  we  get  absorption  at  the  boundary  and  no  transmission  through 
material. 


^  Landau,  L.,  and  E.  F.  Lifshitz.  Electrodynamics  of  Continuous  Media  (New  York:  Pergamon,  1960). 

®  Veselago,  V.  G.  January-Febmary  1968.  The  electrodynamics  of  substances  with  simultaneously 
negative  values  ofe  and  p,.  Soviet  Physics  Uspekhi  10  (4):  509-514.  Available  at 
http://www.iop.Org/EJ/article/0038-5670/10/4/R04/PFlU  10  4  R04.pdf?request-id=be94154d-cd29- 

4796-aca6-209e7c9b79bb 
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Figure  2-1.  An  EM  Beam  Interacting  With  Four  Possible  Materials 

We  now  consider  the  conditions  that  lead  to  negative  permittivity  (e)  and  negative 
permeability  (//)  in  materials.  Negative  values  occur  when  resonances  for  the  electric  and 
the  magnetic  processes  are  in  the  materials.  Negative  refraction  occurs  when  these  reson¬ 
ances  are  at  or  near  the  same  frequencies. 

In  most  naturally  occurring  materials  and  over  most  frequencies,  s  >  0  and 
ju  >  0.  However,  when  resonances  occur  in  a  material,  both  s  <  0  and  ju  <  0  can  occur 
over  a  small  range  of  frequencies.  Generally,  natural  materials  with  e  <  0  occur  over  a 
region  of  the  spectrum  higher  in  frequency  than  those  with  ju  <  0.  Thus,  the  frequency 
ranges  are  far  apart,  and  they  do  not  overlap,  as  shown  in  Figure  2-2.  They  are  not 
expected  to  overlap  in  one  material. 

Negative  s  occurs  in  metals  in  the  optical  region  and  in  semiconductors  in  the 
terahertz  to  infrared  (IR)  region.  Negative  jx  occurs  at  much  lower  frequencies,  tailing  off 
toward  the  terahertz  and  IR  regions.  Examples  are  as  follows: 

•  e  negative  in  gold,  silver,  aluminum,  and  so  forth 

•  jx  negative  in  ferromagnetics  or  antiferromagnetics. 


2-5 


jL/<<?  £<0 

Figure  2-3.  EM  Spectrum  Showing  Frequency  Windows 
for  Potential  Location  of  Negative  ij  and  £  in  Natural  Materials 

Natural  NIMs  (e  and  fL  both  negative)  are  rare,  and  early  researehers  thought  they 
did  not  exist  at  all.  However,  some  were  reported  recently.  Pimenov  et.  al.’^  stated,  “It  is 
generally  believed  that  nature  does  not  provide  materials  with  negative  refraction  ...”  but 
then  the  authors  went  on  to  describe  just  such  a  material:  colossal  magnetoresistance 
manganite  (La2/3Cai/3Mn03)  in  a  high  magnetic  field.  The  effect  was  demonstrated  in 
2007. 

Another  naturally  occurring  NIM  is  bismuth,  which  is  active  at  a  wavelength  of 
2  =  60  pm.8  One  wonders  how  many  others  are  to  be  discovered.  These  materials  may 
have  a  high  attenuation;  therefore,  they  may  not  be  useful  for  transmission.  However,  the 
index  of  refraction  is  still  negative. 


^  Pimenov,  A.,  A.  Loidl,  K.  Gehrke,  V.  Moshnyaga,  and  K.  Samwer.  “Negative  Refraction  Observed  in 
a  Metallic  Ferromagnet  in  the  Gigahertz  Frequency  Range.”  Phys.  Rev  Lett.  98  (19)  (2007):  197401-1- 
197401-4. 

^  Podolskiy,  V.  A.,  L.  V.  Alekseev,  and  E.  E.  Narimanov.  “Strongly  Anisotropic  Media:  The  TFlz 
Perspectives  of  Left-Flanded  Materials.”  Jourwa/  of  Modem  Optics  52  (16)  (2005):  2343-2349. 


2-6 


3.  Basic  Physics  Phenomena 


3,1  Radiating  Dipole 

In  modeling  the  interaetion  of  DNM,  we  are  eoneerned  with  two  proeesses: 

1.  How  to  model  the  absorption  and  emission  of  individual  atoms,  moleeules,  or 
other  entities 

2.  How  to  model  the  transport  of  EM  radiation  through  a  real  material. 

The  purpose  of  this  modeling  is  not  to  provide  answers  to  eomplex  real  situations 
but,  instead,  to  show  the  phenomena  in  its  simplest  form.  For  the  interaetion  of  individual 
atoms  with  the  radiation  field,  we  will  assume  an  eleetrie  dipole  interaeting  with  the  EM 
field  (see  Figure  3-1)  and  model  it  as  a  foreed  harmonie  oseillation  with  damping.  The 
transport  of  radiation  through  a  material  will  assume  a  ehange  in  veloeity  that  results 
from  a  ehange  in  the  density  (i.e.,  index  of  refraetion). 

A  radiating  dipole  produees  the  pattern  of  an  equal  plane  of  radiation,  as  shown  in 
Figure  3-1.  An  EM  field  induees  an  absorption  by  the  dipole.  This  absorption  is  followed 
by  a  reradiation  of  the  EM  energy.  If  the  material  is  eomposed  of  atoms  that  have  dipoles, 
the  radiation  is  absorbed  by  the  dipoles,  whieh  subsequently  reradiate  (elastie  proeess)  or 
absorb  (inelastie  proeess)  and  store  the  energy.  These  phenomena  ean  be  modeled  as 
foreed  oseillations  of  a  mass  on  a  spring  with  drag. 


Figure  3-1.  Electric  Dipole  Radiation  Pattern 
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3,2,  Resonances  in  Materials  (Forced  Oscillations  -  Mass  on  a  Spring) 

Resonance  phenomena  are  exhibited  by  a  simple  model  of  a  charge  on  a  spring 
that  is  moving  under  the  influence  of  an  oscillating  electric  field.  Figure  3-2  shows  a 
small  mass  attached  to  an  immovable  wall.  For  our  demonstration,  this  depiction  models 
an  electron  attached  to  a  much  heavier  atom  under  the  influence  of  an  oscillating  force — 
an  EM  wave. 


Figure  3-2.  Mass  on  a  Spring  Attached  to  a  Wall 


Consider  the  motion  of  the  small  mass,  m,  influenced  by  an  oscillating  force 
Fcos(ccft)  governed  by 


where: 


(Or, 


( d  \ 

m — YX{t)  +  b  — X{t)  +  kX{t)  =  FQCOs{cot) 
dt  \dt 


b 


(3.1) 


—  ,  y  = - ,  and  F^  =  -eE  represent  the  resonant  frequency,  the  damping 


m  2m 

term,  and  the  magnitude  of  the  force,  respectively. 


The  solution  for  the  position  is 


X  = 


where, /((Z»,(Z»o)  = 


V  m  j 
(c0o"-co^y+4y^co^ 


f{(0,  (Uo)cos((ut)  + 


I  m  ) 


and  g{(o,(OQ)  =  f 


g(®,®o)sin(®0> 


(3.2) 


{(o^^-a^f+Ay^o)^ 


and,  as  we  shall  see,  we  can  define  the  absorption  amplitude  as  Aab  = 

V  m  j 


I  m  , 


g((i),(Og)  and 


the  elastic  amplitude  as  Ae/  = 


f((0,(0o)- 


The  amplitude  of  the  elastic  term  Ae/  =  (-eE/m)f  (co,coo)  and  the  absorptive  term 
Aab  =(-eE/m)g  (co,coo)  is  shown  graphically  in  Figure  3-3.  The  constant  Aab  is  called  the 
absorptive  amplitude,  and  the  constant  Aei  is  called  the  elastic  (or  dispersive)  amplitude 
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Figure  3-3.  Absorptive  (Biue)  and  Dispersive  (Green) 

Ampiitudes  (Reai  and  imaginary  Components  in  Some  Deveiopments 

Using  Compiex  Notation)  for  the  Response  From  a  Forced  Harmonic  System 

because  of  the  way  they  deal  with  the  input  energy.  The  time-averaged  input  power 
absorption  is  entirely  due  to  the  Aahsin(cot)  term.  The  Aeicos(a)t)  term  contributes  to  the 
instantaneous  power  absorption,  P(t),  but  averages  to  zero  over  one  cycle  of  steady-state 
oscillation. 

Therefore,  since  this  energy  is  not  absorbed,  it  must  be  transmitted  or  elastically 
scattered.  These  results  follow  from  the  fact  that  the  instantaneous  power,  P(t),  is  the 
force  Focos(cot)  times  the  velocity  dx(t)/dt  =  (coAahCos(cot)  -  coAeisin(cot)).  The  instantane¬ 
ous  velocity  has  a  contribution  that  is  in  phase  with  the  force  and  a  contribution  that  is 
90  deg  out  of  phase  with  the  force.  Only  the  velocity  contribution  that  is  in  phase  with  the 
force  contributes  to  the  time-averaged  power,  P.  This  “in-phase  velocity”  is  contributed 
by  the  “out-of-phase”  displacement,  Aahsin(cot). 

This  result  is  summarized  as  follows: 


Force: 

F{t)  =  Fq  cos(a>t) 

(3.3) 

Position-distance: 

X{t)  =  A^f^  sin(iMf)  -h  A^j  cosifot) 

(3.4) 

Velocity: 

=  -coA^i,  cos{cot)  +  coA^i  sin(n;f) 
dt 

(3.5) 
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•  Instantaneous  power  absorbed  (which  is  given  by  the  product  of  force  times 
velocity): 

P{t)  =  F(t)(dX  /  dt)  =  FQCcA^i^(cos(cot)  cos(cot))  -  FQcaA^i(cos(at)sin(cot))  (3.6) 

•  Average  power  over  one  cycle: 


2  1 
P  -  FqCoA^^  <  cos  o)t  >  -FgCoA^,  <  cos  ft;? sin  cot  >  -  —Ff^coA^^  -  0 . 


(3.7) 


The  medium  is  said  to  have  normal  dispersion  when  resonances  are  present  in  the 
material  (e.g.,  at  a  frequency  coq)  and  when  to  the  left  of  the  resonance  (co  <  coq),  the 
index  of  refraction  is  increasing.  The  medium  is  said  to  have  anomalous  dispersion 
(which  has  been  observed  in  atmospheric  studies)  when  resonances  are  present  in  the 
material  (e.g.,  at  a  frequency  coq)  and  when  to  the  right  of  the  resonance  (co>  coo),  the 
index  of  refraction  is  decreasing.  In  some  materials,  the  index  of  refraction  could  even 
become  negative,  which  brings  us  to  the  topic  of  this  study.  We  begin  by  looking  first  at 
resonance  phenomena  and  then  at  the  derivation  of  the  index  of  refraction  and  its  effect 
on  EM  waves  transported  through  a  medium  composed  of  such  resonances. 


3.3  Calculation  of  the  Index  of  Refraction:  Transport  Through  a  Plate 

Let  us  assume  a  source  at  point  S  in  vacuum  that  produces  an  oscillating  EM  field 
at  point  S  given  by 

E,=E„e'*> 

and  at  point  P  is  given  by 

(3.8) 

Inserting  a  plate  of  optically  dense  material  of  thickness  Az  changes  the  field  at  P  by 
slowing  the  velocity  though  the  plate.  The  process  of  transport  through  the  plate  involves 
absorption  and  reemission  of  the  EM  field  by  resonances  in  the  plate.  One  such  event  is 
pictured  in  the  material  in  the  plate  shown  Eigure  3-4,  but  many  events  in  the  path  of  the 
beam  contribute  to  the  process. 

As  shown  in  Eigure  3-4,  the  source,  S,  induces  oscillation  in  the  dipole  in  the 
material  at  point  z’.  Emissions  from  the  oscillating  dipole  at  point  z’  occur  to  the  right 
and  to  the  left  from  the  dipole.  However,  we  are  only  concerned  with  the  emission  to  the 
right  because  this  emission  contributes  to  the  field  at  point  P,  which  changes  from  the 


3-4 


Figure  3-4.  Transport  of  EM  Radiation  Through  a  Material  Showing  the  Source  (S),  a  Dipole 
Absorber/Emitter  in  the  Material,  and  an  External  Point  (P)  at  Which  the  Field  Is  Recorded 


field  in  Eq.  (3.8)  to  that  shown  in  Eq.  (3.9)  because  of  the  time  delay  At  (Eq.  (3.10)) 
caused  by  the  dense  plate  of  width  Az.  The  total  field  at  point  P  is  given  by  the  sum  of  the 
source  term  (Eq.  (3.8))  and  the  field  produced  by  secondary  sources  in  the  plate: 


E 


P 


_  p  ico{t  -  At  -  z  I  c) 


(3.9) 


where 


At  =  {n  -l)Az/c , 


(3.10) 


which  comes  from  replacing  the  travel  time  through  the  region  without  the  plate  Az  /  c 
with  the  travel  time  through  the  plate  and  nAz  /  c  lengthened  by  the  index  of  refraction  of 
the  plate,  n. 

Expanding  Eq.  (3.9)  using  e’’^  ~  1-x  +  .  .  .  etc.  gives 

-  ioin-\){Azlc)E/“^‘^^'^\  (3.11) 

(source  term)  (secondary  sources  in  plate), 


which  is  the  source  field  and  the  field  from  the  secondary  source,  clearly  exhibited  in 
Eq.  (3.11). 


We  now  go  back  to  the  analysis  performed  in  Section  3.2  using  real  functions  and 
reexpress  with  complex  functions  for  which  the  analysis  is  more  transparent.  Again,  the 
position  of  the  charge  is 


X{t)  = 


e-Ea 


m{oi>l  -Q)  +  iy  •  co) 


(3.12) 
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and  the  velocity  of  the  charge  is 


dX{t)  . 

- =  z  •  ry - 2 - 2 - • 

dt  m(a)Q-a)  +iy-ci>) 


(3.13) 


Using  the  result  that  a  plane  of  oscillating  charges  with  charge  density  rj  (charges 
per  unit  area),  after  integration  over  the  plane  produces  a  field  at  point  z  given  by  (see 
Feynman’s  Volume  I)^ 


— 

Ej.p(t)  = - [velocity of  charges at(t-z/c)]  (3.14) 

IsqC 


and  inserting  that  field  from  Eqs.  (3.12)  and  (3.13),  we  get 

E„(t)  =  ^(eE„lm)-  ,  ‘'f  ,  (3.15) 

ISqC  {cOq-co  +iyco) 


Now  comparing  Eq.  (3.11)  with  Eq.  (3.15)  and  realizing  that  they  represent  the 
same  field,  we  get 

(n-l)Az  = - - •  (3-16) 

Is^micoQ  -  CO  +iyco) 


After  solving  for  n,  realizing  that  NAz  =r/  and  including  the  sum  over  possible 
different  resonances  k  so  that  Nk  is  the  volume  density  of  k  resonance,  we  get  for  the 
complex  index  of  refraction 


n 


1  + 


IsqM 


■z 

k 


CO  , 


Nk 

CO  ^  i  7 


(3.17) 


Sorting  out  the  index  into  real  and  imaginary  parts, 

n  =  n  -in  , 


(3.18) 


we  get  an  expression  for  the  transmitted  EM  field  that  clearly  reveals  the  absorption  and 
elastic  scattering  parts; 


El=e-o^"A^/Ce 

[ absorption] 


-ico{n 


-l)Az/Cp  ico{t-zlc) 
[ elastic  scattering] 


(3.19) 


9 


Feynman,  R.  P.,  R.  B.  Leighton,  and  M.  Sands.  Quantum  physics.  Vol.  1  of  The  Feynman  Lectures  on 
(Addison- Wesley,  1964). 
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The  field  at  point  P,  after  passing  through  the  plate,  is  given  by  the  product  of  a 
decaying  exponential  indicating  absorption  and  a  term  representing  the  original  source  in 
front  of  the  plate  and  the  fields  from  secondary  sources  in  the  plate  caused  by  elastic 
scattering.  These  solutions  imply  that  negative  values  of  n  are  possible  in  theory  and  are 
found  at  frequencies  near  resonances.  Also,  even  though  we  have  been  talking  about  res¬ 
ponses  of  atoms,  we  can  replace  the  “atoms”  in  this  system  with  man-made  structures 
that  give  similar  resonant  effects,  such  as  the  split-ring  resonator  (SRR),  nanowires,  and 
so  forth.  Arrays  of  these  structures  will  still  have  a  total  macroscopic  effect  of  changing 
the  velocity  of  EM  radiation.  These  structures  are  examined  in  Section  5. 
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4.  New  Phenomena  and  Applications — 
or  What  Can  We  Do  Now  That  We  Could  Not  Do  Before? 


Department  of  Defense  (DoD)  technology  developers  and  users  are  always 
looking  for  more  effective  and  lighter  components  and  devices  for  all  conceivable  appli¬ 
cations.  Metamaterials,  or  artificially  engineered  composite  structures,  have  recently  been 
receiving  increasing  attention.  These  novel  artificial  materials  have  numerous  potential 
applications  in  science,  technology,  and  medicine — especially  if  their  novel  behavior  can 
be  extended  to  the  technologically  critical  near-fR  and  visible  regions. 

Engineers  hope  to  exploit  the  extraordinary  EM  response  that  DNM  offer.  DNM 
exhibit  a  unique  ability  to  bend  and  focus  light  in  ways  that  no  conventional  materials 
can.  The  great  potential  for  enabling  several  innovative  structures  will  be  explored  in  this 
section.  To  study  the  effect  of  these  materials  on  EM  propagation  in  various  applications, 
we  developed  a  computer  model.  Appendix  D  gives  the  theoretical  basis  for  this  model, 
and  Appendix  E  shows  some  results  obtained  with  the  model. 

The  broad  range  commercial  and  DoD  applications  of  these  materials  is  contin¬ 
gent  on  extending  the  frequency  range  from  microwaves  to  the  optical  region,  broadening 
the  frequency  range  of  operability,  and  improving  the  understanding  of  the  physics  of 
their  EM  transport  properties,  all  of  which  should  lead  to  significant  enhancements  in  the 
properties  of  these  materials  (bandwidth,  loss,  operational  frequency,  and  so  forth).  How¬ 
ever,  we  should  mention  that,  at  present,  there  are  major  limitations  with  many  meta¬ 
material  designs.  Eor  example; 

•  Aspect  angle  dependence  (not  three-dimensional  (3D)) 

•  Polarization-dependent  behavior 

•  Eossy 

•  Narrow  band. 

This  section  is  divided  into  two  parts 

1 .  Phenomenology  of  the  DNM  or  the  new  and  unexpected  physics 

2.  interesting  new  things  we  can  build  with  DNM. 
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Here,  we  diseuss  some  applieations  that  people  thought  of,  others  that  people  just 
dreamed  about,  and  still  others  that  were  too  strange  to  mention  until  now.  Appendix  F 
presents  a  more  eomplete  list  of  praetieal  applieations. 

4,1  Reverse  Cerenkov  Radiation 

An  effeet  expeeted  in  double  negative  (DNG)  materials  is  reversed  Cerenkov 
radiation.  Cerenkov  radiation  in  normal  materials  is  the  result  of  eharged  partieles  trav¬ 
eling  faster  than  the  phase  veloeity  of  light  in  that  material. 

Figure  4-1  shows  wave  fronts  (eireles  that  represent  points  having  the  same 
phase)  from  moving  eharged  partieles  when  they  travel  slower  than  the  speed  of  light  in  a 
normal  material  (a)  and  when  they  travel  faster  than  the  speed  of  light  in  the  normal 
material  (b)  (eireles  represent  points  that  have  the  same  phase).  The  polarization  of  the 
moleeules  in  the  eorresponding  two  oases  is  shown  in  Figure  4-2(a)  and  4-2(b),  where  in 
normal  media  (a),  the  polarization  is  more  symmetrio  around  the  oharge  and  in  (b)  when 
the  oharge  is  moving  faster  than  the  speed  of  light  in  the  medium,  the  polarization  falls 
behind  the  oharge. 

Similar  to  the  sonio  boom  observed  at  supersonio  speeds,  a  ooherent  wave-front  of 
light  is  emitted  at  a  oharaoteristio  angle  given  by  cos(9)  =  l/(li  n),  where  P  =  v/c  and  n  is 


Figure  4-1.  Wave  Fronts  for  Charged  Particles  Moving 
(a)  Slower  and  (b)  Faster  Than  the  Speed  of  Light  in  a  Normal  Material 


If  we  oonsider  a  eharged  partiole  traveling  in  a  DNG  material,  n  is  negative  and 
the  oharaoteristio  angle  of  the  wave-front  is  obtuse.  In  other  words,  the  Cerenkov 
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Figure  4-2.  Polarization  of  Molecules  in  Normal 
Media  for  the  Two  Cases  Shown  in  Figure  4-1  (a)  and  4-1  (b) 
Note  for  Figure  4-2:  Arrows  show  the  velocity  direction  of  the  particle. 


Figure  4-3.  Geometry  of  Normal  Cerenkov  Radiation 


radiation  is  emitted  in  the  “baekward”  direetion.  Several  analyses  on  Cerenkov  radiation 
prediet  higher  signals  in  DNG  materials  than  in  normal  materials.  Many  speeulate  using 
this  effeet  to  ereate  better  Cerenkov  deteetors  with  DNG  materials. 

Figure  4-4  eompares  Cerenkov  wave  fronts  in  normal  refraetion  materials  (a)  and 
DNM  (b).  In  Figure  4-4(a),  since  the  phase  is  traveling  in  the  same  direction  as  the 
energy,  the  wavelets  from  different  positions  combine  so  the  wave  front  formed  in  the 
Cerenkov  process  is  moving  in  the  forward  direction.  In  Figure  4-4  (b),  the  wavelets  from 
different  positions  combine  so  the  wave  front  is  moving  in  the  backward  direction. 
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Figure  4-4.  (a)  Normal  Cerenkov  Wave  and 
(b)  Reverse  Cerenkov  Wave  Due  To  Transport  Through  a  DNM 


4,2  Reverse  Doppler  Effect 

A  remarkable  eonsequence  of  negative  phase  veloeities  in  DNG  materials  is  a 
reverse  Doppler  effect  (see  Figure  4-5).  In  normal,  right-handed  media  (a),  the  phase 
velocity  is  parallel  to  the  group  velocity.  Therefore,  when  an  observer  approaches  a 
radiating  source,  the  observed  frequency  is  higher  than  the  emitted  frequency.  However, 
in  left-handed  media  (b),  the  phase  velocity  is  in  the  opposite  direction.  Therefore, 
instead  of  the  “bunching”  effect  of  wave  fronts,  an  observer  approaching  a  radiating 
source  observes  a  lower  frequency  than  the  emitted  frequency.  This  effect  was  predicted 
by  Veselago  and  has  been  verified  in  several  simulations  and  experiments  involving  left- 
handed  transmission  linesio  and  photonic  crystals.  However,  the  engineering  of  a  reverse 
Doppler  effect  experiment  in  DNG  materials  still  proves  to  be  difficult. 


^  incident  wave  co, 

Awl 

reflected  wave  m  <  co.  Receding 
Boundary 

1  incident  wave  m, 

v/vw| 

mj\f\ 

—  1-^ 

reflected  wave  o  ><b,  Receding 

Boundary 

(a) 

(b) 

Figure  4-5.  Doppler  Effect  on  Receding  Boundary 


Lai,  Anthony,  Christophe  Caloz,  and  Tatsuo  Itoh.  “Composite  Right/Left-Handed  Transmission  Line 
Metamaterials.”  IEEE  Microwave  Magazine  5  (3)  (September  2004):  34-50. 
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4.3  Reverse  Casimir  Effect 


The  Casimir  Effect  is  a  force  between  two  conductive  surfaces  in  a  vacuum 
caused  by  the  difference  in  density  of  vacuum  fluctuations!!  between  the  inside  surfaces 
and  outside  the  surfaces  (see  Figure  4-6).  This  difference  reflects  the  fact  that  the  fre¬ 
quency  of  fluctuations  between  the  surfaces  is  limited  to  those  virtual  photons  whose 
wavelengths  fit  a  whole  number  of  times  into  the  gap.  Thus,  the  effect  is  to  bring  the 
plates  together  so  the  force  is  attractive. 


Figure  4-6.  Difference  in  Vacuum  Fluctuations 
Between  Plates  and  an  External  Region  for  Normal  Materials 


The  attractive  Casimir  force,  Fcasimir,  between  two  plates  of  area  A  separated  by  a 
distance  L  has  been  calculated; 


Casimir 


TT'h'C  . 
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480 


(4.1) 


where  h  is  Planck's  constant  and  c  is  the  speed  of  light. 


We  see  from  this  equation  that  the  plates  have  to  be  close  together  to  get  an 


appreciable  effect.  Assuming  L  =  1  pm  and  A  =  1  m^ 


F, 


3.14 -6.6 •10^"''  -S-lOAi 


Casimir 
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(4.2) 


The  gravitational  force  or  weight  of  1  kg  is  about  10  N,  so  1.3  mN  is  the  weight  of 
0.13  grams — small  but  measurable.  Plates  that  have  an  area  of  1  cm  at  the  same 


! !  A  concept  in  quantum  field  theory  where  the  virtual  particles,  photons,  appear  and  disappear  out  of  and 
into  vacuum.  Phenomena  such  as  spontaneous  emission,  the  Lamb  Shift,  and  the  Casimir  Effect,  can 
be  explained  in  this  way) 
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separation  would  feel  a  force  equivalent  to  the  weight  of  about  10'^  grams  or  a  force  of 
approximately  lO'^N.  This  force  is  quite  small  but  is  now  within  the  range  of  modern 
laboratory  force  measurement  techniques.  This  force  has  been  measured  by  various 
research  groups,  and  even  detailed  measurements  investigating  the  effect  as  a  function  of 
the  thickness  and  material  properties  of  the  plates  have  been  accomplished.  12 

On  the  other  hand,  inserting  a  DNM  between  the  plates  produces  a  reverse  effect, 
and  the  force  becomes  repulsive,  as  suggested  by  several  theoretical  investigations.! 3, 1 4 
The  negative-index  metamaterial  (according  to  researchers)  is  able  to  modify  the  zero- 
point  fluctuations  between  the  surfaces  by  focusing  the  photons  and  reversing  the  direc¬ 
tion  of  the  Casimir  force.  This  repulsive  force  is  strong  enough  (according  to  researchers) 
to  levitate  an  aluminum  mirror  500-nm  thick,  causing  it  to  hover  above  a  perfect  lens 
placed  over  a  conducting  plate.  At  the  present  time,  we  do  not  have  any  experimental 
verification  of  this  reverse  effect.  Much  controversy  and  much  interest  have  been  gener¬ 
ated  in  the  community.  This  concept  is  theoretical — no  experimental  observations  yet. 

This  reverse  Casimir  effect  could  be  very  important  in  the  fabrication  and  appli¬ 
cation  of  small  systems  because,  as  small  as  it  is,  the  normal  Casimir  effect  contributes  to 
viscosity  in  nanosystems  and  causes  parts  in  nanosystems  and  Microelectromechanical 
Systems  (MEMS)  to  stick  together.  The  reverse  effect  could  be  used  to  alleviate  the 
problems. 

4,4  Reverse  Snell’s  Law 

In  a  normal  transparent  (dielectric)  material,  an  EM  beam  is  refracted  to  the  right 
as  shown  in  Eigure  1-2.  The  angle  of  refraction  is  obtained  from  Snell’s  law,  which  states 
that  the  ratio  of  the  sines  of  the  angles  of  incidence  and  refraction  is  equal  to  the  ratio  of 
velocities  in  the  two  media  or,  equivalently,  to  the  inverse  ratio  of  the  indices  of  refrac¬ 
tion  ni  and  n2.  Thus, 

nisin(6i)  =  n2sin(62),  (4.3) 


!2  Capasso,  F.,  J.  N.  Munday,  D.  lannuzzi,  and  Ft.  B.  Chan.  “Casimir  Forces  and  Quantum 
Electrodynamical  Torques:  Physics  and  Nanomechanics.”  Selected  Topics  in  Quantum  Electr.,  IEEE 
Journal  of  U  (2)  (2007):  400-414. 

!3  Leonhardt,  U.,  and  T.  G.  Philbin.  “Quantum  Levitation  by  Left-Ftanded  Metamaterials.”  New  J.  Phys. 
9  (254)  (2007):  1-11. 

Kenneth,  O.,  I.  Klich,  A.  Mann,  and  M  Revzen.  “Repulsive  Casimir  Forces.”  Phys.  Rev.  Lett.  89  (3) 
(15  July  2002):  033001-1-033001-4. 
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where  6i  and  62  are  the  angles  of  ineidenee  and  refraetion,  respeetively. 

In  a  DNM,  the  law  still  holds,  but,  since  n2  <  0,  and  62  <  0,  n2sin(92)  is  still  posi¬ 
tive,  the  beam  direction  is  to  the  left,  as  shown  in  Figure  1-2. 

4,5  Metamaterials  as  a  Perfect  Lens 

Probably  the  most  talked  about  and  controversial  application  of  DNM  the  ability 
to  focus  incident  waves  from  a  source  beyond  the  diffraction  limit.  In  other  words,  DNM 
can  be  used  to  create  perfect  lenses,  with  no  theoretical  limit  on  resolution. 

The  idea  of  a  perfect  lens  was  first  proposed  by  Pendry.i^  To  summarize  Pendry’s 
argument,  consider  a  radiating  source  in  a  vacuum  and  a  lens  with  an  axis  parallel  to  the 
z-axis.  For  a  given  emitted  frequency,  co,  the  emitted  radiation  will  have  the  following 
Fourier  composition; 

^  Ea{k^ ,  ky  )  exp(z(k,x  +  ky  y  +  k^z  -  cot)) ,  (4.4) 

where  k  =  kxX  +  kvy  +  h,  z  is  the  wave  vector  of  the  Fourier  component,  and  a  denotes 
the  polarization.  The  wave  vector  components  and  the  frequency,  co,  satisfy  the  free- 
space  (positive  refractive  index)  dispersion  relation.  This  constrains  the  value  of  ky 

k^=^co^c^-kl-kl.  (4.5) 

According  to  Eq.  (4.5),  the  maximum  real  value  for  kz  is  \k\  =  coc,  where  c  is  the 
speed  of  light  in  the  medium.  However,  Maxwell’s  equations  allow  for  solutions  where 
kx^  +  ky  >  (coc/.  In  these  cases,  kz  is  positive  imaginary  and  it  is  possible  for  \kz\  >  \k\. 
These  modes  of  emission  are  known  as  evanescent  waves.  From  Eq.  (4.4),  positive 
imaginary  kz  implies  the  amplitude  of  evanescent  modes  decay  exponentially.  Therefore, 
to  examine  modes  of  emission  for  \kz\  >  \k\  (i.e.,  subwavelength  modes),  the  lens  needs  to 
be  very  close  (i.e.,  subwavelength  distances)  to  the  source.  Also,  for  normal  lenses,  these 
modes  will  still  decay  in  the  lens  once  transmitted.  Thus,  evanescent  modes  are  essen¬ 
tially  ignored  and  the  theoretical  limit  on  the  wavelength  of  radiation  that  can  be  focused 
by  a  normal  lens  is  determined  by  kmax  =  \k\: 


deWolf,  D.  A.  “Transmission  of  Evanescent  Wave  Modes  Through  a  Slab  of  DNG  Material.”  IEEE 
Transactions  on  Antennas  and  Propagation  53  (1)  (January  2005):  270274.  Available  at 
http://ieeexplore.ieee.org/ stamp/ stamp.  isp?amumber=0 1377601 
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(4.6) 


.  l-n  l-ivc 
Aw - = - 

^  max  (O 

Now  suppose  the  lens  is  a  DNM  metamaterial  that  is  negatively  matched  to 
vacuum  (i.e.,  =  yUr=  “!)•  According  to  Pendry,  the  information  from  the  evanescent 

waves  can,  in  theory,  be  recovered  with  such  a  lens.  As  discussed  in  previous  sections, 
DNM  are  left-handed  systems,  where  the  dispersion  relation  gives  a  phase  velocity  (or 
the  wave  vector)  that  is  anti-parallel  to  the  group  velocity.  Therefore,  if  the  group  veloc¬ 
ity  does  not  change  at  the  interface  of  the  lens,  the  wave  vector  inside  the  lens  is  the  neg¬ 
ative  of  the  wave  vector  in  vacuum,  (i.e.,  k=  -  k).  This  sign  reversal  also  applies  to 
evanescent  modes.  With  kz  =  —kz,  the  evanescent  modes  exponentially  grow  instead  of 
decay  inside  the  lens,  which  implies  that  a  DNM  can  amplify  evanescent  modes,  and, 
therefore,  subwavelength  information  can  be  recovered  from  the  source,  as  shown  in  Fig¬ 
ure  4-7. 
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Figure  4-7.  Simple  Illustration  of  Perfect  Lensing 

Note  for  Figure  4-7:  The  top  half  represents  a  normal  lens  where  evanescent  waves  continue  to 
decay.  The  bottom  half  represents  a  DNM  lens,  where  evanescent  waves  are  amplified. 

The  amplification  of  the  evanescent  wave  in  the  DNM  (as  shown  in  Figure  4-7) 
has  been  discussed  in  the  literature  by  people  who  argue  that  this  amplification  leads  to 
unphysical  results.  Flowever,  it  also  has  been  defended  on  the  basis  (1)  that  natural  atten¬ 
uation  in  real  materials  presents  the  divergences  in  ideal  situations  from  occurring  or 
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(2)  that  when  the  time  development  of  the  wave  in  the  material  is  examined,  divergenees 
are  suppressed, 

Aeeording  to  Pendry’s  argument,  a  perfeet  lens  with  arbitrarily  small  resolution 
ean  be  eonstrueted  with  DNM,  and  applieations  of  sueh  perfeet  lenses  would  be  vast. 
However,  several  aeademies  are  skeptieal  of  these  elaims.  The  most  noted  erities  in  the 
literature  are  P.  Valanju  at  University  of  Texas-Austin  and  N.  Gareia  and  M.  Nieto- 
Vesperinas  of  Spain's  Higher  Couneil  for  Seientifie  Researeh  in  Madrid.  Valanju  ehal- 
lenges  not  only  Pendry,  but  also  Veselago’s  original  paper  on  negative  refraetion. 
Valanju  argues  that  negative  refraetion  violates  eausality  and  requires  information  to  tra¬ 
vel  faster  than  the  speed  of  light.  Gareia  and  Nieto-Vesperinas,  however,  are  in  agree¬ 
ment  with  Veselago’s  results  but  argue  that  a  perfeet  lens  would  require  infinite  energy 
density. 

Many  experiments  have  been  motivated  by  this  debate,  and  several  authors  have 
reported  observing  perfeet  lens  effeets  in  DNM  for  small  frequeney  intervals.  An  exam¬ 
ple  of  sueh  an  observation  was  reported  by  Lagarkov  and  Kissel, who  report  being 
able  to  image  separately  two  radiating  radio  sourees  that  are  separated  by  less  than  a 
wavelength  using  a  DNM. 

Despite  the  supporting  experimental  results,  many  are  still  skeptieal  of  the  idea  of 
a  perfeet  lens.  This  skeptieism  is  mainly  due  to  the  many  hurdles  physieists  and  engineers 
must  overeome  to  eonstruet  a  praetieal  working  perfeet  lens.  Currently,  the  main  limita¬ 
tion  in  all  reported  perfeet  lenses  is  the  limited  frequeney  range  of  effeet.  Also,  an  ideal, 
perfeet  lens  would  have  isotropie  DNM  properties,  whieh  many  researehers  are  finding 
diffieult  to  eonstruet.  Attenuation  of  signals  is  also  a  problem  with  eurrent  DNM. 


Gomez-santos,  G.  “Universal  Features  of  the  Time  Evolution  of  Evaneseent  Modes  in  a  Left-Flanded 
Perfeet  Lens.”  Phys.  Rev.  Lett.  90  (7)  (21  February  2003):  077401-1-077401-4.  Available  at 
http://prola.aps.org/pdf/PRL/v90/i7/e077401 

Lagarkov,  A.  N.,  and  V.  N.  Kissel.  “Near-Perfeet  Imaging  in  a  Foeusing  System  Based  on  a  Left- 
Handed-Material  Plate.”  Phys.  Rev.  Lett.  92  (7)  (2004):  077401-1-077401-4. 

Kissell,  V.  N.,  and  A.  N.  Lagarkov.  “Superresolution  in  Left-Handed  Composite  Struetures:  From 
Homogenization  to  a  Detailed  Eleetrodynamic  Deseription.”  Phys.  Rev.  B.  72  (8)  (2005):  085111-1- 
085111-8. 
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4,6  Cloaking/Invisibility 


A  potential  application  of  metamaterials  that  has  recently  gained  a  great  deal  of 
attention  is  the  ability  to  create  EM  cloaking  materials.  According  to  Pendry  et  al.,i9  if 
given  a  metamaterial  with  flexible  EM  properties  (i.e.,  control  over  the  permittivity  and 
permeability  tensors),  light  can  be  forced  to  wrap  around  an  enclosed  space  without  scat¬ 
tering  (see  Eigure  4-8).  In  other  words,  the  enclosed  space  is  rendered  invisible.  Several 
simulations  have  supported  Pendry’s  claim,  and  Schurig  et  al.^o  have  achieved  two- 
dimensional  (2D)  cloaking  in  the  microwave  regime  for  a  copper  cylinder. 


Figure  4-8.  An  Artist’s  Rendition  of  EM  Waves  Traveling  Around  an  Enclosed  Space 

Metamaterials  that  would  serve  such  a  purpose  are  not  necessarily  DNG  mate¬ 
rials,  but  anisotropic  media  with  very  complex  material  properties.  Eabrication  of  such 
complex  metamaterials  would  take  advantage  of  the  many  experimental  advances  made 
in  fabricating  DNG  materials,  so  these  subjects  are  related.  The  emerging  field  of  tai¬ 
loring  £  and  p  to  suit  an  application  like  cloaking  is  called  “transformation  optics.” 


Pendry,  J.  B.,  D.  Schurig,  and  D.  R  Smith.  “Controlling  Electromagnetic  Fields.”  Science  312  (7) 
(June  2006):  1780-1782. 

Schurig,  D.,  J.  J.  Mock,  B.  J.  Justice,  S.  A.  Cummer,  J.  B.  Pendry,  A.  F.  Starr,  and  D.  R.  Smith. 
“Metamaterial  Electromagnetic  Cloak  at  Microwave  Frequencies.”  Science  314  (5801)  (10  November 
2006):  977-980. 
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5.  Manufacturing  an  NIM 


To  construct  DNM  or  NIR  materials,  we  need  to  have  eleetrie  and  magnetie 
resonanees  oeeur  at  the  same  frequency.  Pendry  et  ah,  suggested  the  use  of  SRRs  for  the 
magnetie  response^i  and  short  wire  anteimas  for  the  eleetrie  response. 22  A  strueture  of 
thin  wires,  as  shown  in  Figure  5-1,  supports  plasmons  with  a  plasma  frequeney  dependent 
only  on  the  geometry  of  the  structure  and  given  by 

(Op=2-7v-c^  l{a^\n{al  r)),  (5.1) 

where  r  is  the  wire  radius.  The  Drude  permittivity  is  given  by 

s{a))  =  s^{\-a)ll  co^).  (5.2) 


Figure  5-1.  Wire  Structure 


For  the  magnetie  response  refer  to  the  strueture  shown  in  Figures  5-2  and  5-3  and 
eq.  5.3.  The  frequeney-dependent  permeability,  11(0),  for  this  structure  is  given  by 


A(ry)  =  1  - 


1  + 


2/cr|  . 

- -i 

corii^ 


3lc 

f 

nco^  In 
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Ih 

V 
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r 


(5.3) 


21  Pendry,  J.  B.,  A.  J.  Holden,  D.  J.  Robbins,  and  W.  J.  Stewart.  “Magnetism  From  Conductors  and 
Enhanced  Nonlinear  Phenomena.” /iiiiii  Trans.  Microwave  Theory  Tech.  47  (11)  (1999):  2075-2081. 

22  Pendry,  J.  B.,  A.  J.  Holden,  W.  J.  Stewart,  and  I.  1.  Youngs.  “Extremely  Low-Frequency  Plasmons  in 
Metallic  Mesostructures.  Phys.  Rev.  Lett.  76  (25)  (1996):  4773-4776. 
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The  parameters  a,  I,  b,  d,  and  r  are  defined  in  Figures  5-2-5-4,  c  is  the  speed  of 
light  in  a  vacuum,  and  ai  is  the  resistance  of  unit  length  measures  around  the  circumfe¬ 
rence  of  the  rings.  The  resonant  frequency,  (Do,  for  this  structure  is  given  by 


(a) 


(b) 


Figure  5-2.  (a)  Plan  View  of  a  Split  Ring  Showing  Definitions  of 
Distances;  (b)  Sequence  of  Split  Rings  Shown  in  Their  Tracking  Sequence 


(5.4) 


Source:  Adapted  from  Pendry,  J.  B.,  A.  J.  Holden,  D.  J.  Robbins,  and  W.  J.  Stewart.  1999.  Mag¬ 
netism  from  conductors  and  enhanced  nonlinear  phenomena.  IEEE  Trans.  Microwave  Theory 
Tech.  47:2075-2084. 


Note  for  Figure  5-2:  Each  split  ring  comprises  two  thin  sheets  of  metal. 


Figure  5-3.  Plan  View  of  a  Split-Ring  Structure  in  a  Square  Array  (Lattice  Spacing  “a”) 

Source:  Pendry,  J.  B.,  A.  J.  Holden,  D.  J.  Robbins,  and  W.  J.  Stewart.  1999.  Magnetism  from 
conductors  and  enhanced  nonlinear  phenomena.  IEEE  Trans.  Microwave  Theory  Tech.  47:2075- 
2084. 
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Figure  5-4.  Constructing  an  isotropic  Materiai 

Source:  Pendry,  J.  B.,  A.  J.  Holden,  D.  J.  Robbins,  and  W.  J.  Stewart.  1999.  Magnetism  from 
conductors  and  enhanced  nonlinear  phenomena.  IEEE  Trans.  Microwave  Theory  Tech.  47:2075- 
2084. 

Note  for  Figure  5-4:  An  isotropic  material  is  constructed  by  repeating  the  unit  cell  on  the  right  in 
three  dimensions. 


Following  Pendry’s  suggestion,  SRRs  and  wire  antennas  were  used  in  early 
attempts  to  produce  DNM  in  the  megahertz  region.  Figure  5-5  shows  one  of  the  early 
models  of  a  DNM  (metamaterial)  used  for  the  megahertz  regime.  Figure  5-6  shows  the 
Boeing  cube,  another  early  structure  used  for  the  gigahertz  regime.  As  Eq.  (5.1)  and 
Eq.  (5.4)  show,  the  resonant  frequency  for  these  structures  increases  as  their  size 
decreases.  To  build  materials  that  exhibit  negative  refraction  at  the  higher  frequencies, 
researchers  eventually  miniaturized  these  structures  to  tens  of  nanometers  in  size.  How¬ 
ever,  manufacturing  double  SRRs  at  this  scale  became  too  difficult  and  costly.  As  a 
result,  single  split  rings  became  more  common.  Other  methods  that  permitted  smaller 
sizes  and  higher  frequencies  were  also  used.  Eventually,  researchers  turned  to  single  split 
rings  and  wires  in  the  nano  regime,  as  shown  in  Eigure  5-7. 


Figure  5-5.  Wire  and  Split-Ring  Structure  Designed  for  the  Megahertz  Range 
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Figure  5-6.  Boeing  Cube  Structure  Designed  for  the  Gigahertz  Range 


Figure  5-7.  Wire  and  Spiit-Ring  Structures  for  Higher  Frequencies 

Eventually  nanorings  and  nanoantennas  reached  their  limit  of  usefulness  because 
the  small  sizes  resulted  in  lower  conductivity.  Even  silver,  with  its  high  electrical  con¬ 
ductivity,  ran  into  problems  at  this  scale.  This  point  is  noteworthy  because  even  though 
SRRs  and  nanowire  antennas  obtain  their  jj.  and  e  characteristics  from  their  geometric 
structures  (see  Eqs.  (5.1)-(5.3)),  high  conductivity  and  low  dissipation  are  implicitly 
required  for  their  utility. 

To  reach  even  higher  frequencies,  other  approaches,  such  as  the  fishnet  structure 
shown  in  Eigure  5-8,  became  favored.  All  these  changing  structures  are  a  trend  to  get  into 
the  optical  regime. 
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Figure  5-8.  Fishnet  Structures  for  Higher  Frequencies 


Figure  5-9  shows  the  trends  in  the  development  of  DNM  to  reach  the  optical 
regime  and  their  experimental  realizations.  Instead  of  split  rings  and  wires,  fishnet 
designs  with  dielectric  plates  between  conductors,  as  shown  in  Figure  5-8  and  to  the  right 
in  Figure  5-9,  became  popular  in  2006.  The  most  recent  theoretical  suggestion  is  to  get 
away  from  rings  and  wires  and  fishnet  designs  and  use  (1)  plasmon  resonant  metal  nano- 
particles,23,24  (2)  nanofabricated  photonic  materials,^^  or  (3)  dielectric  nanoparticles 
coated  with  conductive  layers.26,27  These  densely  packed  particle  designs  are  expected  to 
provide  not  only  DNM  properties  in  the  optical  regime  but  also  produce  3D  behavior. 


23  Klar,  T.  A.,  Kildishev,  A.  V.,  Drachev,  V.  P.,  and  Shalaev,  V.  M.  “Negative  Index  Materials:  Going 
Optieal.”  Selected  Topics  in  Quantum  Electronics,  IEEE  Journal  of  12  (6)  (November-Deeember 
2006):  1106-1115. 

24  Valentine,  Jason,  Shuang  Zhang,  Thomas  Zentgraf,  Erik  Ulin-Avila,  Denteho  A.  Genov,  Guy  Bartal, 
and  Xiang  Zhang.  “Three  Dimensional  Optieal  Metamaterial  With  a  Negative  Refraetive  Index.” 
Nature  455  (18  September  2008):  376-379. 

23  Lezee,  H.  J,  Dionne,  J.  A.,  Atwater,  H.  A.  “Negative  Refraetion  at  Visible  Frequeneies.”  Science  316 
(5823)  (20  April  2007):  430-432. 

26  Yannopapas,  V.  “Artifieial  Magnetism  and  Negative  Refraetive  Index  in  Three-Dimensional 
Metamterials  of  Spherieal  Partieles  at  Near-Infrared  and  Visible  Frequeneies.”  Appl.  Phys.  A  87  (2) 
(2007):  259-264. 

22  Alu,  A.,  and  N.  Engheta.  “Polarizabilities  and  Effeetive  Parameters  for  Colleetions  of  Spherical 
Nanoparticles  Formed  by  Pairs  Of  Concentric  Double-Negative  (DNG),  Single-Negative  (SNG)  and/or 
Double-Positive  (DPS)  Metamaterial  Layers.”  J.  Appl.  Phys  97  (9)  (May  2005):  094310-1-094310-12. 
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Evolution  of  Structures 
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Figure  5-9.  Trends  in  the  Development  of  Metamaterials 
From  Split  Rings  to  Wire-Mesh  Fishnet  to  Dielectric  Nanoparticles 
Source:  Adapted  from  Soukoulis,  C.  M.,  S.  Linden,  and  M.  Wegener.  January  2007.  “Negative 
refractive  index  at  optical  wavelengths  negative  refractive  index  at  optical  wavelengths.”  Science 
315  (5808):  47-49. 

Note  for  Figure  5-9.  The  solid  symbols  denote  n  <  0,  the  open  symbols  denote  p  <  0.  Insets 
show  fabricated  structures  in  different  frequency  regimes.  Color  scheme:  double  SRRs  (orange), 
U-shaped  SRR  (green),  metallic  nanorods  (blue),  fishnet  structures  (red),  nanospheres  (yellow 
coated). 
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Appendix  A. 

October  2008  Briefing:  Double  Negative  Materials  (DNMs) 


This  appendix  is  on  a  compact  disk  (CD)  attached  to  the  inside  back  cover  of  this 
document.  The  briefing  on  DNMs  was  presented  in  October  2008.  It  contains  movie  clips  (pages 
A-7,  A-29,  and  A-30)  that  show  a  simulation  of  EM  waves  propagating  through  normal  material 
and  through  DNM. 


Double  Negative  Materials  (DNMs) 

Bohdan  Balko 
John  Franklin 
John  Biddle 
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Appendix  B. 

Glossary  for  Double  Negative  Materials  (DNMs) 


Abbe  limit  (diffraction  limit):  A  limitation  of  conventional  optics.  The  resolution  (i.e.,  level  of 
detail)  that  sueh  optieal  systems  ean  achieve  is  limited  by  diffraetion  to  approximately  half  the 
wavelength  of  the  light  being  used. 

Babinet’s  principle:  A  law  of  diffraction  stating  that  the  diffraction  pattern  generated  by  an 
opaque  object  is  the  same  as  that  generated  by  a  hole  of  the  same  size  and  shape,  except  for  the 
intensity  of  the  pattern. 

Cerenkov  effect:  An  effect  that  arises  when  a  charged  particle,  such  as  an  electron,  moves 
through  a  medium  at  a  speed  faster  than  the  speed  of  light  in  that  medium.  (Relativity  dictates 
only  that  a  particle  must  travel  with  a  speed  less  than  or  equal  to  that  of  light  in  a  vacuum, 
whereas  in  a  non-vacuum,  light  may  be  slowed  down  more  than  the  particle.  When  a  particle 
moves  through  a  medium  with  faster-than-light  speed,  it  emits  radiation  known  as  Cerenkov 
radiation.)  The  direction  of  this  radiation  is  different  in  negative-index  materials  vs.  positive- 
index  materials. 

Doppler  effect:  A  change  in  the  frequency  of  a  wave  due  to  relative  motion  between  the  source 
and  the  observer  of  the  wave. 

Far  field:  The  region  whose  distance  from  a  source  of  radiation  is  many  times  the  wavelength  of 
that  radiation. 

Left-handed  material:  Conventional  materials,  which  have  a  positive  index  of  refraction,  obey 
a  physical  principle  known  as  the  “right-hand  rule.”  Thus,  these  conventional  materials  are  called 
“right-handed.”  In  contrast,  negative-index  materials  are  called  “left-handed.” 

Metamaterial:  A  material  whose  microscopic  or  nanoscopic  structure  is  designed  to  produce 
unusual  behaviors  at  a  macroscopic  level.  For  example,  a  negative-index  metamaterial  can  be 
constructed  using  electromagnetic  (EM)  resonators  that  are  much  smaller  than  the  wavelength  of 
light  for  which  negative  index  is  desired.  At  the  size  scale  of  the  light  waves,  the  resonators  are 
“invisible”  (i.e.,  only  the  bulk  behavior  of  the  material  can  be  seen). 

Near  field:  The  region  close  to  a  source  of  radiation,  typically  within  one  wavelength  of  the 
source. 

Near-field  superlens  (NFSL):  A  general  superlens  could  be  constructed  from  a  material  that  has 
negative  refractive  index  (i.e.,  negative  electric  permittivity  and  negative  magnetic  permeability). 
A  special  case  is  the  NFSL.  The  principle  underlying  the  NFSL  is  that  in  the  near  field,  negative 
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permittivity  alone  is  sufficient  to  achieve  the  superlensing  effect.  Thus,  an  NFSL  could  be  con¬ 
structed  from  a  negative-permittivity  film  (e.g.,  silver).  This  type  of  superlens  requires  that  the 
film  be  sufficiently  thin  and  that  the  object  is  sufficiently  close  to  the  lens. 

Negative-index  material:  A  material  that  has  a  negative  index  of  refraction  (i.e.,  that  simulta¬ 
neously  has  negative  electric  permittivity  and  negative  magnetic  permeability). 

Non-linear  optics  (NLO):  Optics  that  employs  media  with  nonlinear  permittivity  (i.e.,  where  the 
permittivity,  8,  varies  with  the  intensity  of  the  light  and,  thus,  the  strength  of  the  electric  field). 
Examples  of  NLO  include  second-harmonic  generation  and  optical  parametric  amplification. 

Optical  parametric  amplification:  A  mechanism  for  generating  a  laser  light  with  a  tunable 
wavelength. 

Perfect  lens:  A  lens  that  can  resolve  images  below  the  Abbe  limit  (see  preceding  definition). 
Such  a  lens  can  be  constructed  using  a  negative-index  material. 

Phase:  A  measure  of  the  distance  from  a  given  reference  point  to  the  nearest  peak  of  a  wave. 

Phase  anisotropy:  A  situation  where  light  waves  with  different  polarizations  have  different 
phases  (see  also  “Phase”  and  “Polarization”). 

Phase  shift:  A  difference  in  phase  between  two  waves  of  the  same  frequency.  Equivalently,  a 
spatial  or  temporal  offset  between  these  waves. 

Polarization:  A  light  wave  is  composed  of  oscillating  electric  and  magnetic  fields.  Polarization 
refers  to  the  directions  in  which  these  fields  point. 

Second-harmonic  generation  (SHG):  A  nonlinear  optical  process  in  which  two  photons  inte¬ 
ract  with  a  nonlinear  material  to  combine  and  produce  a  single  photon  of  twice  the  frequency. 

Superlens:  See  “Perfect  lens.” 
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Appendix  D. 

Computer  Model  of  Electromagnetic  (EM) 

Propagation  in  Double  Negative  Materials  (DNMs) 

D,1  Computer  Model 

In  our  algorithm,  we  use  the  material  model  introduced  in  1900  by  Paul  Drude. 
Other  investigators  use  variations  on  modeling  this  phenomena,  but,  for  this  study,  we 
feel  the  “lossy  Drude  model”  is  sufficient.  Drude  assumed  that  the  interaction  of  a 
material  with  EM  waves  was  well  modeled  by  a  charge  on  a  spring  attached  to  a  wall  or  a 
much  heavier  mass,  subject  to  a  damping  force  and  to  an  oscillating  force.  Resonances 
arising  from  this  interaction  produce  the  interesting  phenomena  that  we  want  to  investi¬ 
gate.  With  these  assumptions,  the  permeability  and  permittivity  are 

(  co^  ^ 

s{a))  =  SQ  1 - 7 — — — r  (U  ,  (U  “electric/magnetic  plasma  frequency”  (D.l) 

and 

co^  '' 

— pm — ^  “electric/magnetic  collision  frequency.”  (D.2) 

a)[a)  +  iYjj 

The  frequency  domain  equations^  for  the  complex  permittivity  and  permeability 
permit  a  frequency  domain  description  of  electric  and  magnetic  current  density: 

2 

P(ry)  =  icos^ - (D.3) 

icoT^  -  CO 

and 

2 

M(n;)  =  ico/u^ - (D.4) 

icoT^  -  CO 


1  Jackson,  J.  D.  Classical  Electrodynamics.  3"^^  ed.  (New  York:  John  Wiley  and  Sons,  1998). 
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One  obtains  the  time-domain  versions  of  these  equations  by  first  rearranging  the 
equations  like  and  then  exploiting  algebraie  properties  of  the  Fourier  transform.  Given 

x(i)  and  X(a))  where  Xfry)  =  i —  f  x(t)exp(-  io)t)dt ,  we  derive  the  Fourier  transform, 

V2;r 

X'(ct)) ,  of  d^x{t)m  terms  of  x(a>) : 

1  1  °° 

X'{a))  =  \im—  -X[o})+  —  \  x{u)exp[-ioM  +  ia)s)du 


■■  lim— [-  x[o))  +  Qxp{ia)s)x{a)y\  =  x((y)lim— [exp(/(yf )- 1] 

^->0  g  £—>0  g 


■■  x{co)\\m— 


g 


E 


{ico)"  i 


(D.5) 


1  00  /  •  \n  n 

■  ^  =  x{co)\\ 


im 

£■^0 


10)  + 


I 


{icoy  s"  ' 


n\ 


■  ICO. 


)X(o) 


Using  the  algebraic  properties  of  the  Fourier  transform,  we  get  from  Eqs.  (D.3) 
and  (D.4)  an  equation  for  Py  and  equations  for  Mnx  and  M^. 


(D.6) 


and 


=  — M/‘normalized  magnetization,”  (D.7) 
Mo 


5«^«z  +  “normalized  magnetization.”  (D.8) 

Mo 

Defining  new  variables,  we  get  the  first  order  system 
Kx  =  d,M^ 

^z  = 

Jy  =  ^>Py 


Combining  the  time-domain  expressions  for  the  currents  expressed  in  terms  of 
these  new  variables  with  Maxwell’s  equations  gives 

VxE  =  -/yo5,H-M 


V  X  H  =  +  P 


(D.9) 


and  yields  scalar  equations  for  the  reduced,  two-dimensional  (2D)  simulation: 
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d,H,^—{d,E,~KX 


(D.10.1) 


d,H,^-  —  (d,E^+KX 

Mo 


d,K^^Mo0^lH^-r„K^, 


(D.10.2) 


(D.10.3) 


(D.10.4) 


(D.10.5) 


S,Jy^^0(^leEy-^e^y 


We  next  form  the  discretized  versions  of  Eq.  (D.IO)  to  get 

H‘:'  y) = ; (i, y) + — {e;  (-,  y) -  (-, j-i'ij-K:  (y y) , 

Mo\_^z 


KT{ij)=  -  \^-^^?\K%j)  +  5,^oColH%j)  ,  (D.11.2) 

i+-r„  LV  2  ;  j 

I  2  "  J 


H‘:\i,j)=H%j)- 


S,  1 


Mo  L^: 


^r'('.y)=  — I U-^7AK%j)  +  S.t,,0>l,H'.(i,j)  ,  (D.11.4) 

l  +  ^Lr  n  2  ;  J 


(D.10.6) 


(D.11.1) 


-[E%j)-E‘^{i-\,j))+K‘Xi,j)  ,  (D.11.3) 


S,  1 


^0  L^z 


[h[  (/, y  + 1)  -  H[ (/, y))  +  ^ [K  (h  j)  -  K {i  (h  j) 


,  (D.11.5) 


^r'(h7)=  (l-^r Jj;(/,y)  +  ^.^o<^;(h7)  .  (D.11.6) 

1  +  ^r  \-\  ^  )  \ 

\  2  X 
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The  execution  of  the  algorithm  requires  the  following  order  for  optimum  effi¬ 
ciency:  (l)’->(3)’->(6)’->(5)’->(2)’->(4)’. 

Figure  D-1  shows  the  submatrix  of  cells  participating  in  the  calculation. 


WJ 

IJ4 

U 

tflj 

Figure  D-1.  Submatrix  of  Cells  Participating  in  the  Calculation 


D,2  Absorbing  Boundary  Conditions 

Due  to  finite  memory  limitations,  computer  simulations  of  EM  fields  are  neces¬ 
sarily  bounded  in  space.  Simple  conditions  on  the  boundaries  will  result  in  the  reflection 
of  scattered  waves,  and  the  reflected  waves  will  interfere  with  incident  and  scattered 
fields  near  the  area  of  interest  in  the  simulation.  For  most  cases,  this  effect  is  undesired. 
In  the  figures  below  (see  Figures  D-2(a)  and  D-2(b)  on  page  D-5)),  a  line  source  is 
allowed  to  radiate  in  a  bounded  region.  At  the  boundaries,  the  first  derivative  of  the  fields 
is  set  to  zero.  After  suitable  amount  of  time,  the  radiated  fields  have  reflected  off  the 
walls  and  propagating  towards  the  source. 

The  simplest  solution  to  resolve  this  problem,  in  terms  of  complexity,  is  to 
enlarge  the  simulated  space.  If  the  space  is  large  enough,  the  fields  in  the  region  of  inter¬ 
est  can  reach  a  steady-state  solution  before  the  scattered  waves  reach  the  boundary.  How¬ 
ever,  enlarging  the  space  will  increase  computation  time  considerably. 

A  more  clever  solution  is  to  implement  an  Absorbing  Boundary  Condition  (ABC) 
algorithm  at  the  boundaries.  As  the  name  implies,  an  ABC  attempts  to  absorb  incoming 
scattered  waves  at  the  boundaries,  hence,  simulating  a  boundless  space.  The  price  that  is 
paid  in  implementing  an  ABC  is  added  code  complexity  and  increased  computation  time; 
however,  with  the  right  ABC,  the  increase  in  computation  time  is  not  considerable. 
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The  ABC  used  in  our  simulations  is  based  on  Liao  et  al.  “multi-transmitting 
theory. ”2  In  Liao's  ABC,  the  values  of  fields  at  the  boundaries  are  approximated  by  extra¬ 
polating  from  the  fields  ealeulated  inside  the  bounded  spaee.  More  preeisely,  eonsider  the 
function  u(x,  t)  and  its  discrete  form  uL: 

Ul  =  u  (Xmax  —  L  h,  t  -  (L  -  1)  At),  (D-12) 

where  Xmax  is  the  position  at  the  boundary  and  h  is  the  spacing  between  elements.  The 
value  of  the  function  at  the  boundary  uo  =  u(xmax,  t  +  At)  can  be  estimated  by  a  poly¬ 
nomial  extrapolation  of  ul  (L  >  0).  Defining  Am  ul  =  Am-1  ul  —  Am-lui  +  1,  where  Al  ul 
=  ul-  ul+1,  an  extrapolated  value  for  ug  of  order  N  can  be  found  to  be 

uq  =  Ui  +  uj  +  A^  U2  +  .  .  .  +zl^  ^  Uj.  (D.13) 

For  N  =  3  (the  order  used  for  our  calculations),  this  simplifies  to 

Uo  =  3ui  -  3u2  +  U3.  (D-14) 

This  extrapolation  method  has  been  shown  to  be  stable  for  h  =  a  cAt  0.5<a<2.  To 
avoid  the  need  for  interpolation  and  to  satisfy  the  finite-distance  time-domain  (FDTD) 
stability  conditions,  we  set  a  =  2.  Figure  D-2(a)  and  Figure  D-2(b)  show  frames  from 
computer  simulations  without  absorbing  boundaries  (D-2(a))  and  when  the  ABC  algo¬ 
rithm  is  used  to  simulate  absorbing  boundaries  (D-2(b)).  Note  the  appearance  of  interfe¬ 
rences  from  reflections  in  (a)  but  not  in  (b). 


Figure  D-2.  Results  Without  (a)  and  With  (b)  Absorbing  Boundaries 


2  Liao,  Z.  P.,  H.  L.  Wong,  B.  P.  Yang,  and  Y.  F.  Yuan.  “A  Transmitting  Boundary  for  Transient  Wave 
Analyses.”  Scientia  Sinica  (Series  A)  27  (10)  (1984):  1063-1076. 
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Appendix  E. 

Modeling  Double  Negative  Materials  (DNMs) 


E,1  Computer  Code  Algorithm 


We  developed  a  computer  code  to  study  wave  propagation  in  DNMs  numerically, 
using  finite-difference  time-domain  (FDTD)  simulations.  To  describe  the  interaction  of 
the  material  with  electromagnetic  (EM)  waves,  we  used  a  “lossy  Drude  model”  of  the 
medium,  which  is  defined  by  a  frequency  dependence  of  the  permittivity  and  permea¬ 
bility  as  follows: 


s{a>)  =  ^’o  1 

V 


a> 


pe 


co{co  +  ) 


(E.l) 


and 


=  1 

V 


CO 


pm 


a){(o  +  irJ 


(E.2) 


where  and  co^^  are  the  electric  and  magnetic  plasma  frequency  andE^  and  r„  are  the 


electric  and  magnetic  collision  frequency. 


We  use  the  following  definitions  to  derive  scalar  equations  for  a  two-dimensional 
(2D)  simulation: 


Mo 

9 

Mo 

9 

9 

d,E^=^[{d^H^-d^HyjX 


(E.3a) 


(E.3b) 


(E.3c) 

(E.3d) 

(E.3e) 


and 
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d,J^=8,CDlE^-TJ^.,  (E.30 

which  are  discretized  and  used  in  the  computer  algorithm,  as  described  in  detail  in 
Appendix  D. 

E,2  Computer  Code  Results 

In  this  section,  we  show  some  simulations  from  our  studies  of  normal  materials 
(n  >  0)  and  DNM  (n  <  0).  These  are  frames  from  movies  that  run  from  the  time  the 
source  is  turned  on  until  a  pre-selected  time.  See  Appendix  A  for  movie  samples  (p.  A-29 
and  A-30). 

Figure  E-l(a)  shows  the  result  of  a  simulation  of  a  wave  propagating  through  a 
rectangular  slab  with  parallel  sides  made  of  normal  material.  For  comparison.  Fig¬ 
ure  E-l(b),  shows  a  wave  propagating  through  a  rectangular  slab  made  of  DNM.  In  (a), 
the  phase  velocity  is  parallel  to  the  group  velocity  (energy  flow).  In  (b),  the  two  velocities 
are  anti-parallel.  Also,  the  wave  of  the  normal  material  (Figure  E-l(a))  is  spreading  out 
while  the  wave  of  the  DNM  (Figure  E-l(b))  is  focusing. 


Figure  E-1.  Propagation  Through  a  Rectangular 
Slab  (Prism)  Made  of  Normal  Material  (a)  and  DNM  (b) 


In  Figure  E-2(a)  and  E-2(b),  the  same  comparison  is  made  for  a  trapezoidal  slab. 
Here,  however,  the  negative  refraction  is  clearly  exhibited  for  DNM  (b)  as  compared  to 
the  normal  refraction  for  normal  material  (a)  and  a  positive  index  of  refraction. 

The  next  set  of  figures  (Figures  E-3a-E-3h)  show  variation  of  refraction  as  a 
function  of  the  index  of  refraction,  n.  The  graph  on  the  right  of  each  figure  shows  the 
assumed  n — a  vertical  line  in  the  figure  intersecting  the  n  curve. 
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Figure  E-2.  Propagation  Through  a  Trapezoidal 
Slab  (Prism)  Made  of  Normal  Material  (a)  and  DNM  (b) 


Figure  E-3.  Refraction  of  an  EM  Beam 
as  It  Passes  From  Air  Into  Materials  With  Different  Index,  n 

Note  1  for  Figure  E-3:  In  (a),  n  =  -1;  in  (b),  n  =  -0.5 


Figure  E-3.  Refraction  of  an  EM  Beam 
as  It  Passes  From  Air  Into  Materials  With  Different  Index,  n 
Note  2  Figure  E-3:  In  (c),  n  =  -0.25;  in  (d),  n  =  0 
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Figure  E-5.  Refraction  of  an  EM  Beam 
as  It  Passes  From  Air  Into  Materials  With  Different  Index,  n 
Note  3  for  Figure  E-3:  In  (e),  n  =  0. 1;  in  (f),  n  =  0.25 


(g) 


(h) 


Figure  E-3.  Refraction  of  an  EM  Beam 
as  It  Passes  From  Air  Into  Materials  With  Different  Index,  n 
Note  4  for  Figure  E-3:  In  (g),  n  =  0.4;  in  (h),  n  =  0.5 
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Appendix  F. 

Metamaterial  Applications 


F.l  Phase  Compensation 

Since  phase  veloeity  is  anti-parallel  with  the  group  veloeity,  a  double -negative 
(DNG)  metamaterial  ean  be  used  for  phase  compensation.  The  example  given  in  Engheta 
and  Ziolkowskii  considers  a  wave  traveling  through  a  double-positive  (DPS)  (normal) 
material  of  thickness  dl  and  index  of  refraetion  nl.  The  wave  will  aceumulate  a  relative 
phase  differenee  after  traveling  through  the  material  dependent  on  nl  and  dl 

nl*k*dl,  (F.l) 

where  k  is  the  wave  number. 

If  the  wave  subsequently  travels  through  a  DNG  material  of  thickness  d2  and 
index  of  refraction  n2  =-\n2\,  the  total  phase  difference  is 

k*(nl*dl  -\n2\*d2).  (F.2) 

Therefore,  the  phase  difference  can  be  eliminated  by  a  proper  choiee  of  n2  and  d2. 

F,2  Dispersion  Compensation  in  Transmission  Line 

Metamaterials  ean  be  used  to  eompensate  for  the  dispersion  of  signals  in  transmis¬ 
sion  lines.  Essentially,  this  involves  ehoosing  e  and  //  such  that  the  eombination  of  the 
uncompensated  transmission  line  channel  and  the  metamaterial  compensator  gives  an 
effective  index  of  refraction  equal  to  1 . 

F,3  Sub  wavelength  Focusing 

This  is  essentially  the  perfect  lens  effect  proposed  by  Pendry^  and  is  already  dis¬ 
cussed  in  the  report. 


1  Engheta,  N.,  and  R.  W.  Ziolkowski,  eds.  Electromagnetic  Materials:  Physics  and  Engineering  Explora¬ 
tions.  (New  Jersey:  John  Wiley  and  Sons  IEEE  Press,  2006). 

2  Pendry,  J.  B.  “Negative  Refraction  Makes  a  Perfect  Lens.”  Phys.  Rev.  Lett.  85  (18)  (2000):  3966-3969. 
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F,4  Zero  Index  of  Refraction 


Construction  of  metamaterials  that  have  small  index-of-refraetion  values  has  been 
reported.  Some  researehers  claim  that  a  zero-index  material  can  be  used  to  produce 
extremely  narrow  antenna  patterns. 

Another  example,  whieh  is  discussed  in  Engheta  and  Ziolkowski,^  considers  an 
infinite  zero-index  cylinder  with  a  line  eurrent  and  surrounded  by  vacuum.  The  solution 
to  Maxwell’s  equations  gives  a  spatially  static  electric  field  inside  the  eylinder  and  an 
oseillatory  solution  outside.  Even  though  the  electric  field  is  statie  inside,  waves  propa¬ 
gate  from  the  cylinder  into  space. 

F.5  Subwavelength  Cavities  and  Waveguides 

Generally,  the  size  and  geometry  of  a  waveguide  determine  the  modes  (frequen¬ 
cies)  allowed  to  propagate  inside.  Also,  in  normal  waveguides,  the  size  of  the  system 
forees  a  frequency  cutoff  below  which  modes  cannot  propagate.  If  the  waveguide  is  filled 
with  a  composite  DPS/DNG  material,  the  cutoff  frequency  can  be  lowered  or  even 
removed  without  changing  the  overall  size  of  the  waveguide.  Also,  for  eertain  geome¬ 
tries,  waveguides  can  be  made  with  metamaterials  to  have  only  one  allowed  propagating 
mode. 


F.6  Resonant  Enhancement  of  Source  Generated  Fields 

Clever  combinations  of  DNG  or  single-negative  (SNG)  materials  with  DPS  mate¬ 
rials  ean  create  ehannels  that  have  an  effeetive  impedance  close  to  zero.  Such  small 
impedance  would  give  large  enhancement  of  souree  emitted  fields. 

F.7  Efficient,  Electrically  Small  Dipole  Antennas 

For  dipole  antennas  surrounded  by  DPS  materials  (or  free  spaee),  the  radiation 
efficieney  of  the  antenna  decreases  signifieantly  with  size.  However,  according  to 
researchers,  if  the  antenna  is  instead  surrounded  by  a  well-chosen  DNG  shell,  the  effi- 
eieney  can  be  dramatically  increased. 

F,8  Artificial  Magnetic  Conductor  (AMC) 

Metamaterials  ean  be  eonstructed  with  low  permitivities  and  large  permeabilities. 
Such  materials  can  be  eonsidered  as  artificial  magnetic  conductors.  Incident  waves  on 


3  Engheta,  N.,  and  R.  W.  Ziolkowski,  eds.  Electromagnetic  Materials:  Physics  and  Engineering  Explora¬ 
tions.  (New  Jersey:  John  Wiley  and  Sons  IEEE  Press,  2006). 
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perfect  AMCs  would  have  zero-phase  reflection.  AMCs  can  be  used  to  resonantly 
enhance  fields  radiated  from  simple  antennas. 

F,9  Beam-Squinting  Reduction  in  Series  Fed  Antenna  Arrays 

Antenna  arrays  typically  require  each  element  to  be  fed  with  signals  in  phase,  and 
makes  use  of  phase-compensated  transmission  lines.  The  phase  compensation  is  fre¬ 
quency  dependent,  so  a  change  in  operating  frequency  leads  to  phase  mismatching  among 
the  antenna  elements,  resulting  in  squinted  beam  patterns.  Theoretically,  this  hurdle  can 
be  overcome  by  using  negative  index  of  refraction  (NIR)  transmission  lines  for  phase 
compensation.  In  theory,  this  phase  compensation  is  not  frequency  dependent,  and,  there¬ 
fore,  beam  squinting  is  reduced. 

F.IO  Backward  Leaky-Wave  Antenna 

Some  antennas  couple  power  in  small  increments  per  unit  length,  either  contin¬ 
uously  or  discretely,  from  a  traveling  wave  structure  to  free  space.  One  example  is  a 
leaky-wave  antenna  that  radiates  microwave  and  millimeter-wave  energy  from  a  series  of 
metal  strips  mounted  on  a  dielectric  waveguide  that  operates  similarly  to  an  optical  fiber. 
The  spacing  of  the  metal  strips  relative  to  the  guided  wavelength  determines  the  direction 
of  radiation.  Leaky-wave  antennas  are  inexpensive  and  easy  to  fabricate,  which  makes 
them  attractive  for  millimeter-wave  systems.  In  the  backward  leaky  wave  antenna,  the 
wave  propagates  away  from  the  antenna,  but  the  radiation  propagates  away  from  the  sur¬ 
face  at  some  angle.  In  the  forward  leaky  wave  antenna,  the  radiation  is  in  the  forward 
direction.  The  angle  of  radiation  depends  on  the  phase  matching  at  the  surface. 

Applications  of  leaky-wave  metamaterials  include  tunable  antennas,  beam¬ 
forming  structures,  inexpensive  complex  conventional  phased  arrays,  compact-range 
wave  generation  and  measurements,  and  endfire  antennas. 

F,ll  Resonance  Cone  Antennas 

A  planar  wire-grid  network  loaded  with  closely  spaced  orthogonal  capacitors  and 
inductors  is  the  metamaterial  used  in  the  resonance  cone  antenna.  It  is  positioned  over  a 
ground  plane  parallel  to  the  plane.  It  is  exited  by  a  single-frequency  point  source  and 
exhibits  conical  high-field  regions  (hence,  the  name).  When  two  such  metamaterials  are 
interfaced,  the  cones  exhibit  negative  refraction  as  well  as  subwavelength  focusing. 

Applications  might  include  instantaneous  spectrum  analyzers  or  antenna 
multiplexers. 
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F,12  Filters  and  Diplexers  (Microwave) 

The  resonance  features  of  DNG  materials  make  them  suitable  for  a  wide  range  of 
filter  and  frequency  multiplexing  applications.  However,  losses  and  asymmetric  fre¬ 
quency  responses  are  still  hurdles  for  this  type  of  application. 

F.13  Tunable  Transmission  Lines 

Transmission  lines  made  from  metamaterials  are  more  tunable  than  their  typical 
counterparts — which  allows  the  construction  of  transmission  lines  of  greater 
functionality. 
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